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We study the large time dynamics of a macroscopically large quantum systems under a sudden 
quench. We show that, first of all, for a generic system in the thermodynamic limit the Gibbs 
distribution correctly captures the large time dynamics of its global observables. In contrast, for 
an integrable system, the generalized Gibbs ensemble captures its global large time dynamics only 
if the system can be thought of as a number of noninteracting uncorrelated fermionic degrees of 
freedom. The conditions for the generalized Gibbs ensemble to capture the large time dynamics of 
local quantities are likely to be far less restrictive, but this question is not systematically addressed 
here. 
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Recently a problem of an evolution of a quantum sys- 
tem after its Hamiltonian suddenly changed (a problem 
of quantum quench) attracted a lot of attention. A par- 
ticular question which arises in this context is whether 
long time asymptotics of the observables of this system 
can be thought of as obeying a thermal Gibbs-like dis- 
tribution The problem can be formulated in the 
following way. Suppose a system is initially in a quan- 
tum state \ip) which is not an eigenstate of its Hamilto- 
nian (but rather an eigenstate of the Hamiltonian be- 
fore the quench). That state can be decomposed as 
= J2 n Cn l n )> wnere \ n ) are eigenstates of the Hamil- 
tonian after the quench. Then the time evolution of an 
expectation value of an observable can be found directly 
as 



(iP\0(t)\i>) = ^2c*c m e 



.(E„-E„ 



^{n\0\m). (1) 



motion of the system are local (which implies that they 
scale with the volume of the system), then there appears 
to be little doubt that any average of this operator in the 
thermodynamic limit (the limit of an infinite volume at 
fixed density) can be described by a (generalized) Gibbs 
distribution. To see that, one can follow any standard 
arguments regarding the behavior of a small subsystem 
interacting with a large bath whose role is played by the 
rest of the system. More detailed arguments specifically 
for the quenched transverse field Ising model are given 
here [6|-l8j. Yet, suppose we are interested in global ob- 
servables. Then it is not obvious that the averages are 
given by the (generalized) Gibbs distribution. In partic- 
ular, one may want to calculate the so-called diagonal 
entropy Q 



Sd = - Pn In Pn 



(3) 



At long times the oscillating terms in this expression are 
supposed to die out (at least average away if one resorts 
to averaging over time, to avoid dealing with situations 
where oscillations may persist forever) , leading to the 
time independent average 

(V>| 0(t) |V) t largc -> P- H O \n) , p n = |c„| 2 , 
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n 

A question arises whether the averages over arbi- 
trary probabilities p n are equivalent to averaging over 
a Gibbs ensemble, with an appropriately tuned temper- 
ature. Moreover, in some cases the system under con- 
sideration is integrable, that is, it has a large number of 
conserved quantities equal to half of the number of its 
degrees of freedom. Then the question is whether the av- 
erages over the arbitrary p n are equivalent to averaging 
over generalized Gibbs ensemble which includes all the 
integrals of motion in addition to energy in its construc- 
tion [mm. 

Before any general discussion, one must emphasize that 
if O is a local operator (in space) , and if the integrals of 



The question may arise whether this entropy is equal 
to the entropy calculated using the (generalized) Gibbs 
ensemble. 

Here we show averaging a nonintegrable system over 
some arbitrary probability distribution p n in the ther- 
modynamic limit is indeed equivalent to averaging over 
a Gibbs ensemble. However, averaging an integrable sys- 
tem over some arbitrary probability distribution is equiv- 
alent to averaging over generalized Gibbs ensemble only 
if the system can be though of having fermion-like de- 
grees of freedom (which appears to be true for a large 
number of integrable systems), and when these degrees 
of freedom are not correlated. We will see, in particu- 
lar, that quenches in the transverse field Ising model are 
not described globally by the generalized Gibbs ensem- 
ble because its effective degrees of freedom, while being 
fermionic, are correlated (while locally the quenches in 
the transverse field Ising model are indeed described by 
the generalized Gibbs ensemble @|). 

The reason why uncorrelated fermionic degrees of free- 
dom are well described by a Gibbs-like ensemble is very 
simple. A single noninteracting fermionic degree of free- 
dom can be thought of as a two level system, one level 



2 



representing an empty fermionic state and one level rep- 
resenting a filled fermionic state. An arbitrary probabil- 
ity distribution of such a two level system is completely 
fixed by its average particle number. Therefore, if the 
distribution produces correct average particle number, it 
is an exact correct distribution. In particular, a Gibbs 
distribution for such a system is exact as well. This gener- 
alizes to a system consisting of a number of uncorrelated 
fermionic degrees of freedom. Conversely, we will see be- 
low that if a system cannot be thought as splitting into 
uncorrelated fermionic degrees of freedom, it appears un- 
likely that it can be described by its generalized Gibbs 
ensemble. 

In what follows we first go over the derivation of the 
equivalence between the microcanonical and Gibbs en- 
sembles for large closed systems, then discuss the diago- 
nal entropy, the generalized Gibbs ensemble, and finally 
go over the example of a quench in the transverse field 
Ising model whose large time behavior does not reduce 
to the generalized Gibbs ensemble. 

Gibbs ensemble. Consider a large system with energy 
levels E n . Its microcanonical entropy is defined by 



e s m (E )= j2s(E-E n )5E, 



(4) 



where 8E is a fixed energy interval. Let us show that 
this definition of the entropy coincides with the canonical 
entropy in the thermodynamic limit. Indeed, 



,S m (E) 

6E 



dp 

— ( 

2tt 



where F((3) is the canonical free energy 



F((3) = -/T 1 In 



-PE„ 



(5) 



(6) 



Free energy is proportional to the volume of the system, 
and so is E. Therefore, the integral over (3 can be taken 
using the saddle point approximation. The saddle point 
equation reads 



The canonical entropy is given by 



dF 

d0' 



(7) 



(8) 



therefore Eq. © really states that F = E - S/j3. 
The solution of that equation (3(E) is the inverse tem- 
perature which corresponds to a specific energy E. 
Once this solution is known, the solution to the inte- 
gral Eq. |(5j) is, within the saddle point approximation, 
exp \J3 (E — F(/3))], so it gives for the microcanonical en- 
tropy 



S m (E) = P(E-F(f3)) = S({3). 



(9) 



This concludes the proof that these entropies are the 
same (more details can be found in Ref. The main 

criterion of the applicability of the proof is that the sys- 
tem is large so the quantity in the exponential in Eq. ([5]) 
is proportional to its volume. Note that the energy in- 
terval SE drops out in the thermodynamic limit and so 
will be omitted in further discussions. 

Diagonal entropy. Let us show that quite generally, the 
diagonal entropy Eq. ([3]) is equal to the canonical entropy 
defined above. Indeed, by virtue of Eq. 0, Eq. is 
equal to 



Sd = - j dEe s ^ +ln ^ In p(E) 



(10) 



(here p(E n ) — p n ), while the normalization condition for 
Pn is 



dEe S m (E)+ln P (E) = L 



(11) 



Entropy is an extensive quantity which scales with the 
volume. In p(E) typically is also proportional to the vol- 
ume (although this condition may be violated in some 
specially constructed p n ). If so, we can take the integral 
over E in both Eq. (fT0|) and Eq. (fTTj) by the saddle point 
approximation. The saddle point occurs when 



dS m (E) l dp(E) 
dE p dE 



= 0. 



(12) 



Eq. (TTT1) gives, at the value E corresponding to the saddle 
point, 



p(E) = e 



- p S m (E) 



(13) 



Finally, the equation Eq. (fTU)) gives, in the same saddle 
point approximation, when combined with Eq. (j9j), 



S d = S m (E) = S((3). 



(14) 



This concludes the proof the the diagonal entropy and the 
Gibbs entropy are equal in the thermodynamic limit, just 
as observed in Ref. [llj in their study of a time evolution 
of a non-integrable system. Let us reiterate that this 
result is obtained assuming that, (a) the logarithm of 
the probabilities In p n are proportional to the volume of 
the system, which is natural (since the total density of 
states typically grows exponentially with the volume) but 
could be violated in some specific cases depending on the 
initial state |^>), (b) p n are smooth functions of energy, 
which again could be violated in some cases and (c) the 
saddle point condition Eq. (fT2"]l has only one solution, 
which again may be violated for some specially chosen 

Pn- 

Generalized Gibbs ensemble. It is straightforward to 
generalize the above reasoning to the case when there are 
extra conserved quantities, such as the particle number 



3 



N. The definition of the microcanonical entropy, replac- 
ing Eq. (g}, is 



,S m (E,N) 



= J2 S (E~ E n )6(N - N n ) SE, 



(15) 



where N n is the number of particles in a given energy 
level. The rest of the formalism changes accordingly, 
with F becoming the grand canonical free energy, the 
integration in Eq. ([5]) being over /3 and the chemical po- 
tential p, and the integral in Eq. (ITU|) getting replaced by 
the integrals over E and N. The end result is the same, 
that is, the diagonal and the Gibbs entropy are equal to 
each other in the thermodynamic limit. 

However, integrable systems have a large number of 
conserved quantities, equal to half of the number of their 
degrees of freedom. Consider for example a Tonks gas 
which maps into a system of one dimensional noninteract- 
ing fermions. Its conserved quantities are the fermionic 
occupation numbers Ik, where k labels one dimensional 
momenta (Ik are real numbers between and 1 while we 
reserve the notation n k for an integer taking values and 
1). Its microcanonical entropy is given by 



j2 n*( 7 *- n *) 



n»,=0,l k 



ra fc =0,l 



ndpk 
9-tt 



. (16) 



where p k are the generalized chemical potentials. Doing 
the sum over n& gives 



Here 



Q(p) 



ln[l 



(17) 



(18) 



n=0,l 



the generalized free energies. The Legendre transform of 
defines the generalized Gibbs entropy by 



s(i) = - £ (vkh + no**)) ,h = - 



dtt(pk) 
dp k 



(19) 



As before, the probabilities to occupy various states 
in this system are given by p(n) where n stands for a 
collection of rij., each rik for each k takes two values, or 
1. Now suppose p(n) factorizes into a product 



p{n) = JJpfe(n fc ), 



(20) 



that is, each degree of freedom is independent random 
variable in the initial state. The normalization condition 
states that 



p k (0)+p k (l) = l 



(21) 



for each k. Then the diagonal entropy is given by the 
sum of the diagonal entropies 

S d = -X>*(1) lnp fc (l) + (1 - p*(l))ln(l - Pfe(l))] ■ 

k 

. (22) 

One can establish by inspection that this entropy is ex- 
actly equal to the Gibbs entropy defined in Eq. (fT9"]l. if 
one identifies 



h — n k pk(rik) — Pfe(l), 



(23) 



that is, Ik is the average occupation number of the initial 
state. Indeed, we write 



h = - 



dVl(pk) 



1 



(24) 



dp k e~^k + i ' 
This allows us to express pk in terms of Ik and evaluate 

s = - + "(/**)) = 

k 

- {Ik ln/fc + (1 - I k ) ln(l - I k )) , (25) 

k 

which of course exactly coincides with the answer for Sd ■ 
Moreover, the probabilities of observing a state with 
an occupation nk in the Gibbs ensemble are given by 

p k 7 (n)=e Q( -^ +nkn . (26) 

It is easy to see that these probabilities give 

p£(0) = l-J fc , p$(l) = h, (27) 

that is, precisely the right values of the probabilities pk as 
given in Eq. ([23)) . That is, the Gibbs ensemble reproduces 
the distribution p(n) exactly, not just for the diagonal 
entropy, but also for all the averages. 

On the one hand, this appears to be nothing but the 
manifestation of the theorem Eq. (|14|) . On the other 
hand, the conditions under which the theorem Eq. (|14p 
was proven no longer hold true. Indeed, by introducing 
a large number of conserved quantities, whose number 
scales with the system size, we have a situation where 
each individual conserved quantity does not scale with 
the system size and the saddle point approximation can- 
not be valid. Let us look at it in more details. 

Eq. (|T7l) allows us to rewrite the diagonal entropy as 



Srf = _W^ lnpfc(J) , 

7.. •* 



\npk(I)-nI-Q 



while the normalization condition Eq. (|2ip reads 



dpdl ln 



(28) 



(29) 
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for each k. Eqs. (1281) and (|29l) constitute the exact def- 
inition of Sd- There is no reason why saddle point ap- 
proximation is valid for the integrals in these equations, 
since the expressions in the exponential are not large in 
the limit of large system size; in fact, they do not depend 
on the size of the system. Nevertheless, computing these 
integrals via the saddle point approximation we find the 
correct answer Eq. (|22p. the same as would have been 
found had these integrals been computed exactly. At the 
same time, the saddle point approximation gives for the 
value of the integral Eq. (|28p precisely S, the generalized 
Gibbs entropy defined in Eq. ([19]) . This appears to be 
the the origin why the Generalized Gibbs ensemble well 
describes the time evolution of integrable systems. 

To see how this works out, we note that evaluating the 
integral Eq. (|2"51) via a saddle point approximation gives 



This means that 



Pk(I) = 



- P -s(i) 



(30) 



(31) 



where S is the entropy calculated within the generalized 
Gibbs ensemble. Substituting this into Eq. (|28l) . we find 
that Sd, evaluated via a saddle point approximation, co- 
incides with S. 

While it is not clear why the saddle point approxima- 
tion gives exact answers in this case, it appears to be 
crucial that we are describing a system which consists 
of a number of independent degrees of freedom, each of 
them taking only two values. If any of these conditions 
are no longer true, the Gibbs ensemble stops describing 
our system. Indeed, the arguments leading to Eq. (|2"T)l 
already hint that it was crucial that we were dealing with 
a degree of freedom which takes only two values. For such 
a variable, if you know its average I you also know the 
probabilities that it takes these two values, which must 
be the origin why the Gibbs distribution is equivalent to 
an arbitrary distribution in this case. 

For example, consider a situation where the initial 
state has p(rik) which correlate some of the rik jl2j. For 



example, suppose p(rii, 712) does not split into a product 
of two terms, depending on n\ and 712 respectively. Then 
the contribution of these two n to the diagonal entropy 
becomes 



d = ~J (2vf ^p(h,I 2 )x 

\np(Ii ,/ 2 )-A 1 l-fl -M2-f2-f2(A 1 l)-n(^ 2 ) 



(32) 



There is no reason why the saddle point approximation 
when applied to this integral should give the same an- 
swer as if this integral is calculated exactly. And in- 
deed, a straightforward evaluation of the saddle point 
approximation shows that the Gibbs entropy and the di- 
agonal entropy are not equal. On the other hand, If 



pih^h) = px{Ii)pi{Ii)i then the saddle point approxi- 
mation gives the correct answer, as this reduces to the 
example of uncorrelatcd degrees of freedom considered 
previously. 

Moreover, the Gibbs entropy Eq. (|19l) in these exam- 
ple of two correlated degrees of freedom by construction 
depends on two parameters I\ and J2, each being the av- 
erage of m and 72.2 respectively, with respect to p. Yet the 
exact entropy depends on three parameters, the entries 
of the matrix p(ni,ri2) constrained by the normalization 
condition. It is possible to change p so that Sd changes 
yet ii, I 2 , and thus the Gibbs entropy S, as well as Gibbs 
probabilities p defined in Eq. (|2l>|) . do not change. Thus 
we see quite generally that Sd and S cannot be equal to 
each other barring some coincidence. 

We deduce from here that the generalized Gibbs en- 
semble correctly describes those systems which can be 
mapped into independent uncorrelated fermion-like de- 
grees of freedom. Introducing any kinds of correlations 
results in corrections to the generalized Gibbs ensemble. 
Those corrections may be small and may not be immedi- 
ately seen in numerical simulations (and may explain the 
discrepancy between the generalized Gibbs and diagonal 
entropies in Ref. 11|). 

A quench in the transverse held Ising model. For an 
example of a quench in a free fermion-like problem with 
strong correlations, consider a quench in the transverse 
field Ising model, which has been a subject of intense 
studies recently [13|- The model is given by the Hamil- 
tonian H = ~J^2jTj — PY^j T j T j+i- We consider a 
situation where the system described this Hamiltonian is 
in its ground state, and subsequently the parameters 7 
and P are changed abruptly. Under the Jordan- Wigner 



transformation rf 



■ aide 



this 



maps into a free fermion problem with the Bogoliubov- 
de-Gennes-like Hamiltonian 



H 



^ hk ia\ak - a_ fc aL fc J + A fe ( a\.(V_ k + a_ fc a fc 



fc>o 



where 



h k = 2(7 + (3 cos k), A fc = 2(3 sin k. 



(33) 



(34) 



This Hamiltonian is diagonalized by a Bogoliubov trans- 
formation from at to bk (known as Bogoliubov particles). 
The ground state is a vacuum of 6-particles, equivalent 
to the BCS ground state. 

If the parameters of the Hamiltonian suddenly change 
the old ground state becomes a linear superposition of the 
eigenstates of the new Hamiltonian, which is diagonalized 
in terms of new "after the quench" fermions Ck ■ It is easy 
to see that this state can be written as 



|G5 id) = ( 



Uk + Vkc' k c'_ k 



10} , 



(35) 



fe>0 



5 



where 



h° k h k 



A° fc A fc 



((Alf + (hlf) (At 



hi) 



(36) 



1 — u\, and hi 



A° k and /i* 



and Afc are expressed in 
terms of the parameters of the problem before the quench 
and after the quench respectively, according to Eq. (f3"4"|) . 

In other words, there are strong correlations between 
particles of momentum k and —k which appear in pairs 
only. It is easy to check now that its diagonal entropy is 
equal to 



Sd 



E 

fc>0 



, lnul+(l-4) ln(l- U 2 fc )] 



(37) 



At the same time, the Gibbs entropy which can be com- 
puted in a straightforward way with the formalism intro- 
duced above in Eq. ([25)1 is 



S 



[u 2 k ln^ + (l-*4)ln(l-*4)] 



(38) 



fe>0 



that is twice as large. It is not equal to the diagonal 
entropy since it does not take into account correlations 
between the free fermions. 

The discussion so far was directly applicable to models 
which maps into free fermions. However it is well known 
that a wide variety of integrable systems behave as a 
collection of fcrmionic degrees of freedom. Typically one 
needs to specify a set of integers, and the values of the 
so-called rapidities fixing the wave function depend on 
whether a particular integer is present in the set, leading 
to a fermion-like degree of freedom n k specifying whether 
a particular integer k is present 14j. Thus everything 
discussed here seems to apply to quite general integrable 
systems. 
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